library(tidyverse)
library(DeclareDesign)
library(rdss)

sims <- 2000
bootstrap_sims <- 2000

source("code/declarations/declaration_15.1.R")

diagnosands <- declare_diagnosands(
  bias = mean(estimate - estimand),
  rmse = sqrt(mean((estimate - estimand) ^ 2))
)
diagnosis_15.1 <- diagnose_design(declaration_15.1, diagnosands = diagnosands, sims = sims, bootstrap_sims = bootstrap_sims) 

write_rds(diagnosis_15.1, file = "diagnosis_objects/diagnosis_15.1.rds")
